Nucleon Excited States in N/=2 lattice QCD 
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We investigate the excited states of the nucleon using Nf — 2 twisted mass gauge configurations 
for pion masses in the range of about 270 MeV to 450 MeV and one ensemble of Nf — 2 Clover 
fermions at almost physical pion mass. We use two different sets of variational bases and study 
the resulting generalized eigenvalue problem. We present results for the two lowest positive and 
negative parity states. 



I. INTRODUCTION 

Understanding the excitation spectrum of hadrons and 
in particular that of the proton is still a challenge. In 
particular, the Pii(1440MeV) positive parity resonance 
known as the Roper, still remains a puzzle having a mass 
lower than the negative parity state S'ii(1535MeV). This 
ordering is contrary to the prediction of constituent quark 
models where the negative parity is lower in mass than 
Pii. Lattice QCD simulations have recently reproduced 
the mass of the low-lying baryon states using gauge con- 
figurations with pions having mass close to the physical 
value [II H] • In these studies volume and cut-off effects 
have been taken into account by performing the calcula- 
tion at different volumes and lattice spacings. Contrary 
to the low-lying baryon states the study of excited states 
has not yet reached the same level of maturity. In or- 
der to extract excited state energies, a robust analysis of 
simulation data keeping systematic errors under control 
is needed. 

The study of excited states is mostly based on the 
variational principle, which was first applied to extract 
glueball masses [3j. One considers a number of inter- 
polating fields as a variational basis and a generalized 
eigenvalue problem (GEVP) is defined, which yields the 
low-lying energy levels. The GEVP has been applied 
recently to study hadron spectroscopy by a number of 
lattice groups [3HH]- A crucial question of such an ap- 
proach is the convergence of the energy levels to the true 
value. This was first addressed in a paper by Liischer and 
Wolff j9j and recently by the Alpha-collaboration [10]. 
In this work, we explore the variational approach as put 
forward by the Alpha-collaboration to study the excited 
states of the nucleon in the positive and negative parity 
channels. We examine two types of nucleon interpolat- 
ing fields as well as different levels of gaussian smear- 
ings. The approach proposed by the Alpha-collaboration 
is compared with the standard GEVP where the refer- 
ence time to is kept fixed at a small value. The main 
outcome of this comparison is that, within the current 
statistical accuracy typically used for baryon calcula- 
tions, namely 0(10^) configurations, we do not see any 
improvements to the standard analysis. Having estab- 
lished at one ensemble of twisted mass fermions that 



the standard approach performs equally well, we adopt 
it for the other ensembles. We analyze a total of five 
ensembles of Nf = 2 twisted mass fermions with pion 
mass in the range of about 270 MeV to 450 MeV and 
lattice spacing a = 0.089 fm determined from the nu- 
cleon mass [2] . Cut-off effects on the mass of the nucleon 
and hyperons were examined in Refs. [21 lllj respectively 
using two smaller lattice spacings. The conclusion was 
that cut-off effects were within the statistical errors and 
one could take the continuum limit assuming negligible 
0{a^) contributions. Therefore, in this work, we limit 
ourselves to studying only one lattice spacing for twisted 
mass fermions. In addition, we analyze an ensemble of 
Nf — 2 Clover fermions with pion mass m^r ^ 160 MeV 
and lattice spacing a ~ 0.073 fm. 

The paper is organized as follows: In section|Tljwe give 
the details of the simulations; in section III we compare 



results using different variational basis and analysis ap- 
proaches for an ensemble of twisted mass fermions with 
TOjr ^ 300 MeV; in section IV we give our results, and in 
section [Vj we summarize our findings and give our con- 
clusions. 



II. SIMULATION DETAILS 

The input parameters of the calculation using Nf = 2 
twisted mass fermions, namely /3, L/a and a/i are sum- 
marized in Table |lj These are the same configurations 
already used in the analysis of the low-lying baryon spec- 
trum |llj , where more details regarding the twisted mass 
formulation can be found. The corresponding lattice 
spacing a and the pion mass values, spanning a mass 
range from 270 MeV to 450 MeV, are taken from Ref. [5]. 
For the lattice spacing we take the value determined from 
the nucleon mass, which, taking into account the system- 
atic error due to the chiral extrapolation, is consistent 
with the one extracted from /^r [T^ . 

Apart from the twisted mass fermion ensembles given 
in Table [T] we also analyzed an ensemble of Nf ~ 2 
Clover fermion configurations produced by the QCDSF 
collaboration. We use the 48"^ x 64 ensemble with near- 
physical pion mass of ~ 160 MeV, a.t f3 — 5.29 
for which the lattice spacing has been determined to 
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P — 3.9, a — 0.089(1)(5) fm from the nucleoli mass 
rp/a ^ 5.22(2) 



24^ X 48, afi 0.0040 0.0064 0.0085 

L = 2.05 fm No. of confs 782 545 348 

(GeV) 0.3131(16) 0.3903(9) 0.4470(12) 

m^L 3.25 4.05 4.63 

32=* X 64, a^j. 0.003 0.004 

L = 2.74 fm No. of confs 659 232 

(GeV) 0.2696(9) 0.3082(6) 

m^L 3.74 4.28 



TABLE I. Input parameters (/3, L, ji) of our lattice calculation 
and corresponding lattice spacing (a) and pion mass (m.^). 
The value of the lattice spacing determined from /,r is a = 
0.0855(6) fm. 



be a = 0.0728(5) (19) fm [IT. This yields a value for 
wLtj L ~ 2.8. We smear the links that enter the Dirac op- 
erator with three iterations of APE smearing and set the 
clover term to its tree-level value i.e. csw = 1- Smear- 
ing the links in this way chang it. We therefore 
tune the value of the hopping parameter k as described 
in |14j to match the pion mass in the unitary theory. 
A comparison of the pion and nucleon effective masses, 
ameff{t) = log [C{t)/C{t + 1)], in the unitary theory and 
after tuning is shown in Fig. [ij As can be seen, the mass 
of the nucleon in the tuned theory agrees with the one ob- 
tained in the unitary theory. Note that one has to allow 
10 time slices or about 0.7 fm to ensure that excited states 
have been sufficiently suppressed. Performing a double 
exponential fit yields a suppression factor of 0(e^^) from 
the gap from the mass gap between the ground and ex- 
cited state, which means there is a substantial overlap 
of the standard nucleon interpolating field with higher 
excited states. 
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FIG. 1. The pion (red circles) and nucleon (blue squares) 
effective masses in the non-unitary setup as described in the 
text, compared to their values in the unitary theory (solid 
black line) computed by QCDSF [T^. The value of k in the 
non-unitary setup was tuned to reproduce the pion mass in 
the unitary theory. 



nucleon interpolating fields, namely 

Xi = {u^Cj5d)u and X2 = {u^Cd)-f^u. (1) 

The local scalar-diquark nucleon interpolator, xii is well 
known to have a good overlap with the ground state of 
the nucleon, while the X2 interpolator, which vanishes in 
the non-relativistic limit, is believed to couple strongly 
to the first excitation of the nucleon, the Roper. In ad- 
dition, the variational basis is enlarged by considering 
different gaussian smearings using similar parameters to 
those used in Ref. |T0] . The correlation matrix considered 
here thus has the general form: 
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III. THE VARIATIONAL METHOD 

The standard extraction of the ground state energy 
from the large time limit of Euclidean two-point correla- 
tion functions relies on the fact that they are expressed 
as a sum of the energy eigenstates of QCD that exponen- 
tially decay as a function of the time with a rate propor- 
tional to the energy. The variational method provides an 
approach for extracting, besides the lowest energy state, 
the low-lying excited states from Euclidean correlation 
functions. A variational basis is constructed by using dif- 
ferent interpolating fields x with the quantum numbers 
of the particular state of interest, which in this work is 
the nucleon. Applying the variational principle one can 
determine the superposition of states that correspond to 
the low-lying nucleon states. One variational basis is ob- 
tained by considering two different spin combinations of 



W - E Trf^ (1 ± 7o) (xi'Hx, Oxi'^ (0, 0))] 

X 

oo ^ 

= ^ e-^"*Tr[- (1 ± 7o) (Olxi^ V) ("Ixi'^ |0)] , 

rt=0 

i,j = l,...,7V 

a,& = 1,2, ^ ' 

where the trace is taken over Dirac indices and C^{t) 
{C~{t)) yields the positive (negative) parity correlator. 
The states \n) are eigenstates of the Hamiltonian with 
En < En+i and we have assumed that the temporal ex- 
tent of the lattice is large enough to neglect boundary 
contributions. We note, however, that in our analysis 
we average over the forward and backward propagating 
nucleons. The indices i and j on the correlation matrix 
[t] correspond to different levels of gaussian smearing 
and a and h to Xi a-nd X2- 
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A. Variational basis with different gaussian 
smearing levels of xi 

In this subsection we perform an analysis using a vari- 
ational basis of xi with a number of different smearing 
levels. The variational basis is constructed using N dif- 
ferent gaussian smearing levels of this interpolating field. 
The GEVP is defined by the generalized eigen-equation 

C{t)Vn{t,to) = Xn{t,to)C{to)Vnit,to), 

n = 1, . . . ,7V, t > to , (3) 

where En — lim(_i.oo — <9f log A„(i, to)- The corrections 
to En decrease exponentially like e""^-^"* where AEn = 
min„j^„ \Em — En\ ^ for fixed Iq. In Ref. 0110] it was 
shown that if one varies to such that to > t/2 then the 
correction is 0(e~^^" "*) with AE^^n = E^ - En en- 
suring a greater rate of convergence. In this section we 
examine the benefit of this relation for extracting the 
low-lying states in the nucleon sector. We apply gaus- 
sian smearing to each quark field, q(x.,t) [13 US], enter- 
ing xi- The smeared quark field is given by q^'^'^^'^{x, t) — 
J^(x, y; U{t))q{y, t) using the gauge invariant smear- 
ing function 



F(x, y; U{t)) = (1 + aHT' (x, y; U{t)), 



(4) 



constructed from the hopping matrix understood as a 
matrix in coordinate, color and spin space 

(x, y; U{t)) = (U,{-K, t)(5x,y-aJ + C// (x-ai, <)(5x,y+a?] . 

i=l ^ ^ 

(5) 

Smearing is applied at the fermion source and sink. 
Following Ref. [TU] we consider values of the smearing pa- 
rameters a = 0.1 and =0, 22, 45, 67 and 135. These 
smearing parameters produce a source with a root mean 
square (r.m.s.) radius in lattice units of 0, 1.96, 2.72, 
3.25 and 4.48 respectively. These different smearing lev- 
els are labeled by the superscript i = 1, ... ,5 on x^*-*- 
The resulting correlation matrices are symmetrized. We 
use 150 twisted mass fermion configurations with /3=3.9, 
a/u = 0.004 or m^ - 308 MeV on a 32^ x 64 lattice. 

Let us first examine the advantage of using these differ- 
ent smearing levels. We consider several different corre- 
lation matrices of the positive parity correlator ^ . (t) 

constructed from Xi for different smearing levels i in 
order to examine both the role of varying and/or the 
dimensionality of the GEVP. In Fig. [2] we show the effec- 
tive mass for the ground and first excited states result- 
ing from a GEVP analysis of all possible 3x3 correlation 
matrices fixing to /a = 1. We are looking for the com- 
bination of interpolating fields that give the fastest con- 
vergence to the two-lowest levels Eo and Ei i.e. to the 
earliest onset of a plateau behavior. From this analysis it 
is evident that using the higher smearing levels slightly 
improves convergence. Nevertheless, the overall quality 
of the plateaus after time-slice i/a = 6 is similar allowing 
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FIG. 2. The effective mass for the ground (-Bo) and first 
excited {E\) states resulting from a 3 x 3 GEVP. A 3 x 3 
correlation matrix was constructed out of different interpo- 
lating fields Xi^ by applying a different number of gaussian 
smearing iterations on xi . The numbers in the legend give the 
combination of the three values of Us used to construct the 
basis. The effective energy levels resulting from a truncated 
3x3 GEVP constructed using Eq. Q are also shown. This 
analysis was carried out using 150 configurations of twisted 
mass fermions at /3=3.9, a/i = 0.004 (m^ ~ 308 MeV) on a 
32^ X 64 lattice. 
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FIG. 3. The effective mass for the ground (n — 0) and first 
(n = 1) excited states from the best choice of 2 x 2 and 
3x3 GEVPs corresponding to the highest level of Gaus- 
sian smearing i.e. from {xi^^Xi^'} f™' the 2x2 GEVP and 
{Xi\x^i\xf^} for the 3x3 GEVP. The test was carried out 
using 150 configurations of the twisted mass ensemble with 
/3=3.9, a^i = 0.004 (m^ ~ 308 MeV) on a 32^ x 64 lattice 



extraction of the mass with all ten possible 3x3 combi- 
nations. Apart from making a choice of the best basis 
by trying different combinations of gaussian smearing we 
also try a truncation scheme where the 5x5 correlation 
matrix is projected to an m x m matrix, C™^'"(t), using 



the m eigenvectors belonging to the m largest eigenvalues 
of C(to) with TO < 5 as follows: 



C 



<NxN 



{to)b ^ Ab, 



NxNi 



k,j 



TO, i ~ I, . . . , N, 



(6) 



where Ajk = Sjke~^^*° is a.n N x N matrix with the 
eigenvalues of C^^^(io) as its diagonal elements and b 
a.n N X N matrix with the corresponding eigenvectors. 
We additionally tried this truncation scheme with various 
values of to/a, namely to/a = l,.--,4 and results obtained 
are statistically equivalent. The resulting effective masses 
extracted from the truncated 3x3 matrix, are plotted in 
Fig. [2] As can be seen, they do not show any improved 
convergence, in particular in comparison with those ex- 
tracted from the optimal 3x3 correlation matrix, i.e. the 
"345" combination, constructed with x'i\x^^ and x'l^- 

The effect of reducing the dimension of the GEVP to 
2x2 can be seen in Fig. |3] The quality of the plateaus for 
the first two states is not affected as compared to those 
extracted using the optimal 3x3 correlation matrix. 

In Fig. ID we compare the GEVP results to the stan- 
dard method of extracting the effective mass using a sin- 
gle interpolating field Xi \ or equivalently to the trivial 
1x1 GEVP in the case of the ground state. As can 

(5) 

be seen, using just the Xi interpolating field yields the 
same quality plateau as that obtained from the 3x3 cor- 
relation matrix analysis. This is to be expected since this 
smearing is optimal for the nucleon ground state. 

For the two lowest states in the positive channel we 
also study the resulting eigenvectors in order to under- 
stand/verify the mixture of the various Xi'^ contributing 
in the best interpolating field for each state. Identify- 

(i) 

ing the optimum combination of Xi extracted from the 
GEVP analysis is useful if one wants to calculate the 
matrix elements of any operator using the optimal in- 
terpolating field that best suppresses the contribution of 
excited states. 

In Fig. [5] the three components Vi , V2 and V3 of the 
eigenvector for the ground and excited states in the pos- 
itive parity channel determined from the optimal 3x3 
correlation matrix constructed from xi with smearing it- 
erations Us — 45,67 and 135 are displayed. The interpo- 
lating field with the maximum overlap with the ground 



state is given by Xoff = ViXi'^ + ^2X1' + ^sXi 
V is the large-time limit of v{to) = limt_>oo ^(i, io)- It 
is evident that most of the contribution comes from the 
components vi and V2, while V3 has very little contribu- 
tion (an order of magnitude less). This is in agreement 
with the previous observation that a 2 x 2 GEVP with 
rig = 67 and 135 is adequate for these first two states. 
The ground state is dominated by the highest smearing 
level (vi) and for the first excited state the contribution 
from V2 increases. 

Let us next vary to as suggested in Ref. [10], where it 
was shown that this lead to an improvement in the deter- 
mination of the ground state by successfully suppressing 
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where 




FIG. 4. The effective mass for the ground state for to/a — 1. 
Results are shown for the GEVP with the best choice of 3x3 
using different gaussian smearings on Xi^ with i = 3, 4, 5 and 
from the correlation function C12 12 and C15 15. 
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FIG. 5. The components of the eigenvector for the ground 
and first excited states at to — 1. Results are shown for a 
GEVP with the optimal 3x3 correlator, Ci;,ij, i,j = 3,4,5. 



excited state contamination for certain mesonic systems. 
In Fig. |6] we show results obtained at fixed to/a = 1 as 
well as results obtained by varying to- Within the sta- 
tistical accuracy of our analysis, we see consistent results 
for the three values of to/a = Ij 3, and 5 considered. Fur- 
thermore, we allow ^0 to vary for every value of t and, in 
particular, we apply the condition to > t/2 as suggested 
in Ref. 10, . We show results for the ground and first 
excited states in the positive parity channel for the case 
to = t/2, where we observe no change in the plateau range 
within the present statistics. For these nucleon states and 
within the present accuracy, this analysis does not show 
an improvement. Our conclusion is that for the low- lying 
nucleon spectrum, where the energy gap is not particu- 
larly small, the variation of that has been shown in 
Ref. ^Oj to reduce the systematic error is not observed 
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here at least within the limitation of our statistics. Thus, 
for our case, one can efficiently extract the energy spec- 
trum from the effective mass plateaus at a fixed value of 
to . In what follows we will &x to/a = 1. 



meff(n) for the ground-state (n = 0) and the first excited 
state (n = 1) and an indication for the second and third 
excited states {n = 2 and n — 3). 
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FIG. 6. (a) The effective mass for the ground state for 
various choices of to- Results are shown for the best choice of 
3x 3 GEVP. (b) The effective mass for the ground and first 
excited states with a fixed value for to (squares) and with the 
condition to = t/2 (circles) for the ground (filled symbols) and 
first excited state (open symbols). Values have been slightly 
shifted in time in order to aid the comparison. 
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FIG. 7. The spectrum when using Xi '' for a 5 x 5 GEVP using 
150 configurations at l3=3.9, a^i = 0.004 (m^ ~ 308MeV) on a 
32'^ X 64 lattice. The basis set of interpolating fields consisted 
of Us — 0, 22, 45, 67 and 135 Gaussian smearing iterations. 
The solid lines and bands show the fitted effective mass and 
jackknife error. 

From the above analysis it is clear that the merit of 
the variational approach lies in the extraction of excited 
states, whereas the ground state is equally well obtained 

(5) 

using the optimally smeared interpolating function Xi ■ 
In Fig. [7] we analyze the 5x5 GEVP to extract the 
nucleon spectrum. Despite the low statistics used for this 
test case, we are able to obtain effective mass plateaus 



B. Combining both xi and X2 

In the preceding subsection we used a variational ba- 
sis constructed from different smearing levels of the xi 
interpolating field. In this section, we extend the in- 
vestigation by combining both Xi' ^ind Xj*^ each with 
two different smearing levels: a local source and ns=135 
with a— 0.1, i.e. using the interpolating fields: x^i \ X2^\ 
Xi and X2 ■ thus have a 4x4 correlation matrix. 
We examine the same ensemble as in the previous sec- 
tion. In Fig. |8] we display a comparison of this 4x4 
basis with the full 5x5 basis of the previous section. 
This is done for the ground and first excited states in the 
positive parity channel. The effective mass plateaus are 
statistically equivalent for both basis sets. In addition we 
compare to the 2x2 problem constructed using the basis 
{x^i \ X^2^}- As can be seen this yields results consistent 
with the 4x4 and 5x5 systems. Furthermore, by ex- 
amining Fig. [9J where we illustrate the effective mass of 
the two lowest states extracted from this 4x4 GEVP to 
those obtained from the correlators Ci^ and 25 i-e. 
using alone x^i^ or X2^ both source and sink, we find 
that the effective masses coincide. For the ground state 

this was already demonstrated in Fig. H The new result 

(5) ■ — ^ ■ 

is that X2 yields the first excited state in the positive 

parity channel. 
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FIG. 8. The effective mass for the ground and first excited 
states for the positive parity channel for /3=3.9, a/^ = 0.004 
on a 32'^ x 64 lattice. The 5x5 system is obtained from C\^i- 
with i,j < 5; 4x4 corresponds to the that obtained using 
Caibj with a, & = 1, 2 and i,j = 1,5 and 2 x 2 to the system 
obtained from Ca^t^ with a,b — 1,2. 150 configuration are 
used for the 5 x 5 GEVP and 280 for the 4 x 4 and the 2x2 



This important conclusion is corroborated by consider- 



ing the results shown in Fig. 10 where we compare the ef 
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FIG. 9. The effective mass for the ground and first excited 
states in the positive channeL The effective masses from the 
interpolators x^i^ ^-nd X^2^ '^^ compared with tlie results ob- 
tained from the same 4x4 GEVP as explained in Fig. [S] 
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FIG. 10. The effective mass for the ground and first excited 
states in the positive channel. The effective masses from the 
interpolators Xi^'' and X2^'' is compared with the results ob- 
tained from the 5x5 GEVP as explained in Fig. [S] 



fective mass extracted from a GEVP analysis of Ci. i (i), 

*j J ^ 5 with the effective mass when using just Xa i.e. 
Can a5 with a = 1,2. It is evident that the first excited 
state obtained from the GEVP is identical to the ground 
state extracted from C25 25 , carrying smaller statistical 
errors, a result that we will use in order to further ex- 
amine the first excited state. This illustrates that X2 
strongly couples to the first excited nucleon state but 
not to the ground state. 

For the case of the negative parity states, it is impor- 
tant to note that the negative parity interpolating opera- 
tor in Eq. [2] has a non-zero overlap with the two particle 
S-wave state which consists of a nucleon and a pion. At 
the physical point, this state has lower energy than the 
negative parity nucleon. To know a priori if and at which 



pion mass the negative parity nucleon and the tt N state 
cross requires knowledge of the pion mass dependence 
of the negative parity nucleon. In Fig. [TT] we show the 
ground and first excited states obtained from a 4 x 4 and 
2x2 GEVP. As in the case of Fig. |8l the 4 x 4 GEVP 
is solved using the basis {x''i\x'i\)<2^\x'2^}j while the 
2x2 using {Xi^\x^2 ^}- As can be seen the ground and 
first excited states are statistically equivalent when ob- 
tained by solving the GEVP of either size. 




FIG. 11. The nucleon ground (filled symbols) and first excited 
state (open symbols) in the negative parity channel, evaluated 
via a 4 X 4 GEVP (black circles) and a 2 x 2 GEVP (blue 
diamonds). The interpolator bases used for these two GEVPs 
are the same as those used in Fig. [8] for the same GEVP size. 
240 configurations were used for this test. 

Having verified that the 2x2 GEVP yields the same 
energies for the ground and first excited states of the neg- 
ative parity as the 4x4 GEVP does, we will hereof use 
the 2x2 GEVP to resolve the ground and first excited 
negative parity states for all other pion masses. Knowing 
which of these is the multi-particle state would require in- 
vestigating the dependence of the two energy levels on the 
lattice volume, which is beyond the scope of this work. 
We therefore compare the two energy states with the sum 
of the nucleon and pion mass, and from this infer which 
is the negative parity nucleon state. Further examples 
of the effective masses obtained for the negative parity 



states using the 2x2 GEVP are given in Figs. 12 and 13 
discussed in the following section. 



IV. THE LOW-LYING NUCLEON SPECTRUM 

In the previous section we have shown that if we are 
interested in the first excited positive parity states of the 
nucleon both the variational analysis and the interpolat- 
ing fields Xi and X2 yield equivalent results for the 
ground and first excited state. Since in the rest of this 
work we will be interested in the first excited state we 
minimize the computational effort and carry out an inves- 
tigation of the ground-state and the first excited state of 
the nucleon using the correlators constructed with xi and 
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X2 with a single smearing level. We note that, to mini- 
mize computational resources, we reuse existing forward 
propagators, namely those used in [TT]. These propaga- 
tors were created with smearing parameters rig =50 and 
a=4.0, which yield an r.m.s. radius consistent with the 
parameters ns=135, a^O.l which were used in the pre- 
ceding section and labeled x'^'*- We nevertheless label 
this new smearing level differently, i.e. as x^^''- Equiva- 
lently, for the negative parity states, we solve the 2x2 
GEVP constructed from C~gj,^ with a,b = 1,2, as already 
discussed. We also note that results presented from here 
on have been obtained with the statistics listed in TableU] 



In Figs. 12 and 13 we show the effective masses for 
both positive and negative parity states, for a twisted 
mass ensemble and for the Clover ensemble analyzed in 
this work. Note that the twisted mass ensemble used 



in Fig. 12 is the same as that used in Fig. 11 however 
with almost three times more statistics; in Fig. [Tl] the 
purpose was to perform a test for GEVPs with different 
bases and therefore less statistics were utilized to speed- 



up the calculation. As can be seen from Figs. 12 and 13 
a plateau region can be identified for both the positive 
and negative parity states and the twisted mass fermion 
and clover ensembles respectively. The effective masses 
for the other twisted mass fermion ensembles display a 
similar behavior. 
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FIG. 12. The effective mass when using the Xi'^ 
interpolating fields for the positive parity states (lower panel, 
= i"^), and the 2 X 2 GEVP using C~^^^ with a,b=l, 2, 

for the negative parity states (top panel, = ^ )■ The 
results are for the ensemble with /J=3.9, ajj, = 0.004 on a 
32^ X 64 lattice. The errors were obtained by the jackknife 
procedure. 

The results for all of the ensembles of Table U and 
the single Clover ensemble are displayed in Fig. [14] For 
the nucleon mass we apply continuum chiral perturba- 
tion theory to extrapolate lattice results to the physical 
pion mass, omitting the Clover point from the fit. We use 
SU(2) heavy baryon chiral perturbation theory (HBxPT) 
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FIG. 13. The effective masses of the ground and first excited 
nucleon states for both positive and negative parity, for the 
Nf = 2 Clover lattice at m„ ~ 160 MeV. The rest of the 
notation is the same as in Fig. [12] 
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FIG. 14. The first positive and negative parity states mea- 
sured in our different gauge ensembles. The twisted mass en- 
sembles are plotted with filled symbols, while the results from 
the single Clover ensemble are denoted with the open sym- 
bols. We show chiral extrapolations for the nucleon ground 
state to the 0{p'^) result of Eq.[7]using only the twisted mass 
fermions data. The dashed line is the resulting curve obtained 
by adding the pion mass to the 0{p^) curve. Physical masses 
for the different states are indicated by the magenta circles. 



to 0{p^) given by 

TO7v('Tl7r) = ITT- 



(0) 
N 



4c^) ml 



(7) 



Since the lattice spacing was fixed using the nucleon mass 
it is no surprise that the curve passes through the phys- 
ical value. However, the fact that it also passes through 
the result obtained with clover fermions is a non-trivial 
confirmation of the pion mass dependence of the nucleon 
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mass. In the figure we also show a curve obtained by 
adding the pion mass to the curve describing the nucleon 
mass. As can be seen, for all pion masses considered 
here, the negative parity ground state is consistent with 
the mass of the pion plus nucleon, indicating that this is 
the two particle tt N state in an S-wave configuration. 

In Figs . [T5| and [T6| we compare the results of this work 
with three other calculations available in the literature. 
Namely, we compare with the result using a Clover im- 
proved fermion action by the CSSM collaboration [19] 
with a ~ 0.09 fm, a calculation using anisotropic Clover 
lattices by the Hadron Spectrum Collaboration [20 with 
spatial lattice spacing = 0.123 fm and a calculation us- 
ing the Chirally Improved Dirac Operator by the Bern- 
Graz-Regensburg (BGR) collaboration [21] and lattice 
spacings between 0.13 and 0.14 fm. We note that the 
lattice spacings for the two latter calculations are no- 
tably larger than those used in this work giving rise to 
issues about cut-off effects. 
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FIG. 15. The positive parity states of this work (filled and 
open squares) compared with results from other groups, that 
include a Nj = 2 -\- 1 Clover improved fermion calculation 
by the CSSM collaboration [T^ (red diamonds), a calcula- 
tion using anisotropic clover lattices by the Hadron Spectrum 
Collaboration [23] (open hexagons) and a calculation using 
the Chirally Improved Dirac Operator by the Berlin-Graz- 
Regensburg (BGR) collaboration (21j (yellow triangles). 

The first observation is that all lattice results are in 
agreement in the case of the mass of the nucleon. Fur- 
thermore, in all cases results using twisted mass and 
Clover fermions are consistent. The second major ob- 
servation is that our data for the first excited state of the 
nucleon in the positive parity channel are too high to be 
identified with the Roper. In fact all lattice data except 
the ones from the CSSM Collaboration are incompati- 
ble with the Roper. Our result using Clover fermions 
at TO^n- ~ 160 MeV remains high and shows no indica- 
tion of converging to the mass of the Roper unlike the 
result from the CSSM collaboration. We do not observe 
a sharp change of slope as we decrease the pion mass as 
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FIG. 16. The negative parity states of this work compared 
with calculations from other groups using the same notation 
as in Fig. [15] 



observed in Ref. [T^ and therefore it is unlikely that this 
discrepancy can be totally due to volume effects, since 
both our Clover calculation and the one in Ref. [THj are 
on lattices of comparable pion mass and physical lattice 
extent. In the negative parity channel our results are 
consistent with the ones from the BGR collaboration. 
We can clearly see that for all pion masses considered 
the negative parity ground state is consistent with air N 
state in an S-wave. To the statistical accuracy available 
to us, the first excited negative parity state appears to 
be converging to A^^(1535), however the errors are too 
large to draw concrete conclusions. It is clear from this 
analysis that extracting the excited states is much more 
challenging and more work is needed to fully understand 
the low-lying spectrum of the nucleon. 



V. CONCLUSIONS 

In this work we have applied the variational method to 
investigate the excited states of the nucleon. Two sets of 
variational bases were used and the analysis of the result- 
ing GEVP was performed using the standard approach 
of fixing as well as by varying <o such that > t/2 
as advocated in Ref. [TU]. Within the current statistical 
accuracy, we found that for the nucleon excited states no 
observable improvement is obtained as compared to fix- 
ing to- Limiting ourselves to the first excited state of the 
nucleon in the positive parity channel suffices to use the 
Xi and X2 interpolating fields with the optimal smearing. 
The former couples strongly to the ground state of the 
nucleon whereas the latter to the first excited positive 
parity state. Using these two interpolating fields we did 
an analysis using Nf ~ 2 dynamical fermion configura- 
tions reaching almost physical pion mass. Even at the 
lightest pion mass of 160 MeV we see no evidence for the 
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Roper. This is in agreement with most other lattice cal- 



culations shown in Fig 15 and may indicate that the in- 
terpolating fields used have small overlap with this state. 
In the negative parity state our result at = 160 MeV 
is consistent with the mass of A^~(1535). 



Addendum: After the completion of this work the 
CSSM collaboration updated the results shown in this 
paper (Figs. 15 and 16 1 in Ref. 
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